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Abstract. Instead of a quantum computer where the fundamental units 
are 2-dimensional qubits, we can consider a quantum computer made up 
of d-dimensional systems. There is a straightforward generalization of the 
class of stabilizer codes to d-dimensional systems, and I will discuss the 
theory of fault-tolerant computation using such codes. I prove that uni- 
versal fault-tolerant computation is possible with any higher-dimensional 
stabilizer code for prime d. 



1 Introduction 

Quantum computation and quantum communications have the potential to ac- 
complish many things that would be difficult or impossible using just classical 
computers and communications. However, quantum data is very vulnerable to 
decoherence and to errors. It is likely that some form of quantum error correc- 
tion will be needed to perform anything beyond the simplest computations with 
a quantum computer. Quantum error-correcting codes |l], g, ||, [|, [| provide one 
of the tools necessary. Such a code can protect quantum data against errors 
occuring during transmission or storage of the data. However, to have a reli- 
able quantum computer, we also need for the computation to be performed in 
a fault-tolerant manner || . Fault-tolerant quantum computation requires a pro- 
tocol that not only maps states of a quantum code to other states of a quantum 
code, but prevents errors from propagating out of control. 

A large group of useful codes was introduced in || and Q. These codes are 
essentially the quantum equivalent of classical linear codes, in that they can be 
easily described and encoded, and that it is easy to measure the error syndrome 
(as in the classical case, it may be difficult to compute the actual error from the 
error syndrome) . The primary complication involved in quantum error correction 
is that it is not only necessary to correct bit flip errors 

x=(y ), a, 



but also phase errors 
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Consequently, it turns out to be useful to look at the group V generated by 
tensor products of these operators. V is called the Pauli group or the extra- 
special group. Stabilizer codes are those codes where the valid codewords are all 
eigenstates of n—k operators in V . These n — k operators generate a 2 n ~ k element 
Abelian group, called the stabilizer S of the code. The set of valid codewords 
forms a 2 fe -dimensional subspace of the full rt-qubit Hilbert space, the coding 
space of the code. The stabilizer is analogous to the parity check matrix of a 
classical linear code. In fact, the set of classical linear binary codes is exactly 
the set of stabilizer codes where everything in the stabilizer is a tensor product 
of Z operators. 

Stabilizer codes are easy to work with because of the structure of the Pauli 
group. 

XZ = -ZX , (3) 

so any two operators in V either commute or they anticommute. If an operator 
E anticommutes with an operator G in the stabilizer S of a code, then E\ip) will 
have eigenvalue —1 for G instead of +1. Therefore, by measuring the eigenvalues 
of the n — k generators of S, we can easily measure the error syndrome, and if 
the code is suitably chosen, identify the error. 

Even if we restrict our attention to stabilizer codes, it is far from obvious 
that we can perform fault-tolerant computation. First of all, we must find some 
operators that map the coding space of the code into itself. Secondly, many of 
these operators will cause errors to spread from one qubit in a block to a different 
qubit in the same block of the code. Therefore, a single error could rapidly grow 
to become many errors within a block, exceeding the code's capability to correct 
them. In order to avoid this, we will further restrict attention to transversal 
operations, that is, operations which only interact qubits from one block with 
corresponding qubits in other blocks. This means that an error occuring in qubit 
number 3 in one block might spread to qubit number 3 in another block, but 
it will never spread back to qubit 2 in the first block. Each block can easily 
correct a single error, so this situation causes no problems. Shor was the first to 
demonstrate a universal set of fault- tolerant gates ||. However, his construction 
only worked for a small class of codes. In J7j], I was able to show that Shor's 
construction could be generalized to any qubit stabilizer code. The proof made 
extensive use of the group of unitary operators that leave the group V invariant 
under conjugation. This group is known as the Clifford group. 

In the classical theory of error-correcting codes, it is often helpful to go 
beyond bits and work with higher-dimensional systems. The same may be true for 
quantum error correction. Instead of a system made up of 2-dimensional qubits, 
we can work with a system composed of <i-dimensional audits. It turns out that 
there is a natural generalization of stabilizer codes to higher-dimensional systems 
|8[ p[ p!o|. In Sec. g, I will present this generalization. Then I will procede to 
generalize the arguments of |7) to show that universal fault-tolerant computation 
is also possible with any of these codes, at least in the case where d is prime. 
Though I will largely focus on prime d, I will also say a little bit about the general 
case along the way. Assume that d is prime unless it is otherwise specified. 



The construction of a universal set of gates used in |]7) consisted of a number 
of steps: 

1. The full Clifford group can be constructed given the controlled-NOT, oper- 
ators in the Pauli group, and measurement of operators in the Pauli group. 

2. For any stabilizer code, we can perform encoded versions of all operators in 
the Pauli group and can measure operators in the Pauli group. 

3. We can perform a CNOT between corresponding encoded qubits in different 
blocks of the code. 

4. We can swap individual encoded qubits from one block of the code to an 
empty block and vice-versa. 

5. We can move an encoded qubit in an otherwise empty block to whatever 
position in the block we desire. At this point, we have established that we 
can perform a CNOT between any pair of encoded qubits in the same or 
different blocks, and therefore can perform the full Clifford group. 

6. Given the full Clifford group, we can perform an additional gate outside the 
Clifford group, such as the Toffoli gate or the tt/8 rotation. This completes 
the universal set of gates. 

Each of the steps in this construction has an analog for higher-dimensional sys- 
tems. However, in this paper, I will present a simplified construction that com- 
bines steps |3] through |5[ This construction gives a direct way to perform the 
generalization of the CNOT between any pair of encoded qudits, whether they 
are in the same or different blocks and whether they are in corresponding or 
different positions within their blocks. This construction can also be used to 
simplify the proof for systems with d = 2, and will likely reduce the required 
overhead for fault-tolerant computation using codes with many qudits per block. 



2 Higher Dimensional Generalization of the Pauli Group 
and Other Structures 

The Pauli group has a natural generalization to higher-dimensional systems.! 
Instead of generating it from the two-dimensional X and Z, we instead generate 
V from tensor products of Xd and Zd, where Xd\j) = \j + 1) and Zd\j) = c^\j), 
where a; is a primitive d-th root of unity. Xd and Zd satisfy the relation 

X d Z d = uJ^ZdXd . (4) 

Both Xd and Zd have order d. The elements of the single-qudit Pauli group have 
the form uj a X r d Z s d , where < r, s < d, and 

{xizd (x\zi) = u st -™ (x\zt) {xizi) . (5) 

2 The group presented in this section is not the only generalization of the Pauli group, 
although it is probably the simplest. See j^] for a more extensive discussion of this 
issue. 



V for n qubits will contain d 2n elements, plus an additional factor ofji for overall 
phase. The elements of V have eigenvalues uj t for r — 0,...,d — la From now 
on, I will suppress the subscript d, and all operations should be taken to be over 
qudits instead of qubits. 

The d-dimensional generalization of the stabilizer S of a code is again just an 
Abelian subgroup of V . The coding space is composed of those states that are 
fixed by all elements of S (when d is even, this actually imposes an additional 
constraint on the overall phase of elements of S). If the stabilizer on n qudits 
has n — k generators, then S will have d n ~ k elements and the coding space will 
consist of k qudits. Note that this last fact need no longer be true when d is not 
prime, and this is the main source of complications in that case. It is unclear 
exactly how to deal with a code that does not encode an integral number of 
qudits. If we stick to codes for which all the generators of the stabilizer have 
order d, the rest of the proof will hold, modulo a question about gates necessary 
to generate the Clifford group. 

If an operator E and M £ S satisfy 

EM = uj a ME , (6) 

then E\ip) will have eigenvalue u> a for M instead of eigenvalue +1, so we can 
detect that error E has occurred by measuring the eigenvalue of M. 

We can see the structure of the coding space by extending the generators of S 
to a complete independent set of commuting operators. When d is not prime, we 
also require these operators to have order d. Such a set will have cardinality n, so 
we can do this by choosing k additional operators Z\, . . . , Z^. These operators 
have the interpretation of the encoded Z operators for the k encoded qudits. We 
can then choose k more operators X\, . . . ,Xk which satisfy the relations 

Xtfj = z^{i^j) (7) 
XiZi = lu 1 ZiXi (8) 
XiM = MX, (VAf e V) . (9) 

The operators Xi then act as the encoded X operators for the k encoded qu- 
dits. The generators of S along with Zi,...,Zk and X\,... ,X/~ then generate 
the group of all Pauli group operators that commute with S. As in the two- 
dimensional case, the operators that commute with S but are not themselves 
in S are precisely the operators that cannot be detected by the quantum code. 
They therefore perform encoded operations on the data. 

The Clifford group is the set of operators that leave V invariant under con- 
jugation. That is, it is the normalizer N(V) of V in the unitary group U(d n ). 
The Clifford group is important for fault-tolerant computation because if we 
perform an operator U on the Hilbert space, the operator UNW has the same 
relationships to states after the transformation as the operator N did before the 
transformation. Therefore, instead of considering transformations of the states 



3 This is true for odd d. For even d, XZ has order 2d, so extra factors of i will be 
necessary, as in the d = 2 case. This aspect is actually simpler for odd d than for d = 2. 



\ip) — > U\ip), we can consider transformations of the operators N — ► UNU' f . 
When U is in the Clifford group and N is in the Pauli group, then UNU^ will 
also be in the Pauli group. Therefore, we can uniquely describe elements of N(V) 
by the permutation they induce on V . The permutation must preserve the group 
structure of V, but is otherwise arbitrary. 

In the two-dimensional group, N(V) was generated by two single-qubit op- 
erators R (the Hadamard transform |j) — > |0) + (— and P (the phase gate 
li) -» ij b')) 5 and thc two-qubit operator CNOT -> + j) mod 2)). In 

d dimensions, i? generalizes to the d-dimensional discrete Fourier transform 

b->-Xy» , (io) 

P generalizes to the d-dimcnsional phase gate 

li}-^°'- 1)/2 li) , (ii) 

and CNOT generalizes to the SUM gate 

mod d) . (12) 

We can describe these operators by their induced transformations on the 
Pauli group. R maps 

X -► Z , (13) 

Z^JiT 1 . (14) 



P maps 



SUM maps 



X->XZ , (15) 
Z^Z . (16) 



, (17) 

J®A"-» I®X , (18) 

Z®I Z®I , (19) 

7 ® Z ^ Z" 1 <S> Z . (20) 

However, it is not clear that these three gates generate the Clifford group. We 
may also need the S gate 

X -> X a , (21) 
Z ^ Z b , (22) 

for all pairs (a, 6), where afr = 1 mod d. On kcts, this gate acts as \j) — > |aj). 
In fact, a single pair (a, 6) is sufficient, as long as a generates the multiplicative 
group ~K* d . I will not give a detailed proof that these gates generate the Clifford 



group, but using the P, R, and S gates, we can get the full one-qudit Clifford 
group. Then a construction similar to that used in |t]] will give the full n-qudit 
Clifford group. The structure is somewhat more complicated when d is not prime, 
and I have not verified that these gates are sufficient for the nonprime case. 

Note that we can fault-tolerantly measure any operator in V that is the ten- 
sor product of Z operators by performing transversal SUM gates from the qudits 
to be measured to an appropriate ancilla. Since we are interested in eigenvalues 
with possible values ui^ for j = 0, . . . , d — 1, the appropriate ancilla state is the 
superposition of all states where the sum of the qudits is mod d. That way, no 
information beyond the eigenvalue of the measured operator will be conveyed. 
We can construct this state by Fourier transforming the state X^jLo \tt '"!))■ 
Following DiVincenzo and Shor we can measure any operator in V by per- 
forming a transversal Clifford group operation C that takes the operator to the 
tensor product of Z's, performing the measurement, and applying C _1 . 

For any stabilizer code, the elements of the fc-qudit encoded Pauli group 
are also elements of the n-qudit unencoded Pauli group, as are the generators 
of the stabilizer. Since we can perform and measure an arbitrary element of 
the unencoded Pauli group, we have shown that for a stabilizer code over d- 
dimensions, we can apply encoded versions of X and Z for all encoded qudits, 
measure the generators of the stabilizer (and therefore perform fault-tolerant 
error correction), and measure all members of the encoded Pauli group fault- 
tolerantly. This provides step || of the proof. 

3 Measurements and Stabilizers 

In 0, it proved very helpful in a number of places to understand how the sta- 
bilizer of a state or subspace changed under measurements. The procedure for 
qubits generalizes easily to higher dimensions. 

First, recall that there is more than one way to choose generators for a given 
stabilizer. Any maximal set of independent operators in the group will suffice. In 
particular, any generator M can be replaced by NM for any N ^ M . Similarly, 
the encoded X and Z operators are only defined up to multiplication by elements 
of S. If we wish to measure an operator A £ V , then the first step is to put the 
stabilizer and X and Z operators in a form so that all the X and Z operators 
commute with A and all but one of the generators of S commutes with A. We can 
do this because if M G S does not commute with A, then M a N will commute 
with A for some a for any N that commutes with M (as do all X's, Z's, and 
generators of S) . We will not need to consider the case where A commutes with 
everything in the stabilizer. 

This is a useful form for the stabilizer because any operator that commutes 
with A is not disturbed by the measurement of A. Therefore, we only need to 
change M when A is measured. Since A £ V, the possible measurement results 
are ui a for a = 0, . . . , d— 1. These result uj a corresponds to applying the projection 
operator 

P a = (l + u- a A + Lj- 2a A 2 + ■■■+ u- {d - 1)a A d -^ jd (23) 



to the state. Assume now that 

MA = ujAM (24) 
(note that when d is prime, this will always be true for some power of M). Then 

d-1 d-1 

MP a M ] = 2^2 u~ ja MA j M^ = ~ d Y J v~ j{a ~ 1) A j = P a -i . (25) 

3=0 j=0 

Thus, if we measure A and get the result a, by applying M a we can produce the 
same state we would have gotten if we had gotten the result 0. 1 will assume below 
that any measurement is followed by such a correction. Once this correction is 
performed, the new state is a +1 eigenvector of A, so A should be added to the 
stabilizer. It is not an eigenvector of M, so M is removed from the stabilizer. 
All of the other generators (which have been put in a form where they commute 
with A) are unchanged. 



4 Gates Derived from SUM 

Suppose we have the ability to perform the SUM gate between any pair of qudits 
in our computer, as well as the ability to perform the Pauli group and to measure 
operators in the Pauli group. I will now show that we can apply the full Clifford 
group to the computer. 

Suppose we consider a single unknown qudit and prepare a second ancilla 
qudit in the state |0). This two-qudit system can be described by the stabilizer 
I <g> Z. The logical Pauli group is generated by X — X <g> / and Z = Z ® I . Now 
perform a SUM gate from the first qudit to the second qudit. The stabilizer is 
Z- 1 <8> Z, X = X <g> X, and Z = Z <g> I. 

Suppose we were now to measure the operator A = I ® XZ . Then M = 
Z~ l ® Z G S and MA = luAM. Therefore, this measurement results in the 
stabilizer I ® XZ, and X = XZ^ 1 <8> XZ and ~Z = Z ® I. We can discard the 
second qudit, and the effective transformation on the first qudit is 

X -> XZ- 1 , (26) 
Z^Z . (27) 

This is the gate P^ 1 . d — 1 iterations of it will produce the P gate. 

Alternatively, we could have prepared the ancilla qudit so that the stabilizer 
of the system began as J ® X, then performed the SUM gate from the second 
qudit to the first. The stabilizer would then be X <S> X, and X = X <g> I and 
Z = Z ® Z^ 1 . Then we measure A = I® XZ^ 1 and choose M = X ® X so that 
MA = ujAM. The final stabilizer is / ® XZ -1 , so we discard the second qudit, 
leaving 



X^X , 
Z -> XZ 



(28) 
(29) 



Call this gate Q. Then R- 1 = XQP~ l Q, and R = R~ 3 . 

Now suppose we again prepare the ancilla in the +1 cigcnstate of X, but 
now perform s SUM gates from the second qudit to the first instead of one. The 
stabilizer is X s ® X, X = X ® /, and Z = Z ® Z~ s . This time we measure 
A = Z <E> I. This results in stabilizer Z <E> I, 

X = (X <g> I) {X s <g X) - ^ 1 = Jig X" 8 " 1 , (30) 

and Z = Z®Z~ S . Therefore, discarding the first qudit leaves the transformation 

X -> X-"" 1 , (31) 
Z -> Z~ s . (32) 

By choosing an appropriate s, we can therefore perform an arbitrary S gate. 
Note that in this case, the data ends up in what was originally the ancilla qudit. 

I have shown how to produce the P, R, and S gates from the SUM gate. 
Therefore, given the SUM gate, we can produce the full Clifford group. This 
completes step [l] of the proof. 

5 Producing the SUM gate for any stabilizer code 

To see how to construct the SUM gate between any pair of encoded qudits, first 
consider two unencoded qudits. Introduce a third qudit in the state |0). The 
stabilizer at this point is I<8> I®Z. Assume we can do Pauli group measurements, 
even entangled ones, and perform operators in the Pauli group. Let us first 
measure the operator I ® X ® X^ 1 . This becomes the stabilizer. The logical 
Pauli group generators are 

Xi = X <g> I g> J , 
X 2 = I <g> X ® J , 
Zi = Z ® I <g> J , 
Z 2 = I <E> Z <g> Z . 

Now measure Z ® 1 ® Z . It becomes the new stabilizer, and 

~X X = X ® x ® x- 1 , 
X 2 = I <8> X ® I , 
~Zx = Z ® I ® I , 
~Z 2 = I ® Z (8) Z . 

Finally, measure I <g> I ® X and discard the last qudit. This leaves us with 

x l= x®x , 
X 2 = I ® x , 
~Z X = Z%1 , 
~z 2 = z^ 1 <g> z 



(33) 
(34) 
(35) 
(36) 



(37) 
(38) 
(39) 
(40) 



(41) 
(42) 
(43) 
(44) 



This we recognize as the transformation induced by the SUM gate, so this series 
of entangled measurements has performed the SUM gate between these two 
qudits. 

Now, to apply this to a quantum code, we just need to be able to measure 
entangled logical Pauli group operators between any pair of encoded qudits. If 
the qudits are in the same block, this is straightforward. For instance, if they 
are in slots i and j, the encoded version of X ® X is just XiXj. This is in the 
Pauli group too, so we know how to measure it. 

If the qudits are in different blocks, it is not much harder. Instead of using 
an a-qudit ancilla state, we use an (a + 6)-qudit ancilla state (where a and b are 
the weights of the operators X j and X j ) , which is again in the superposition of 
all states whose registers sum to mod d. The operator we wish to measure is 
Xi®Xj, which is in the Pauli group. By performing the appropriate transversal 
Clifford group operation, we rotate this to be the tensor product of Z's and 
perform SUM gates from the appropriate qudits to the corresponding ancilla 
qudits, then perform the inverse Clifford group operator to restore the state to 
its original form. Then we measure the a + b ancilla qudits, and this tells us the 
eigenvalue of the measured operator. We use an (a + 6)-qudit ancilla instead of 
an a-qudit ancilla plus a 6-qudit ancilla because we do not wish to be able to 
find the eigenvalues of Xi ® I and / Cg> Xj separately, only their product. 

Therefore, given an encoded ancilla qudit which is initialized to |0), by per- 
forming the encoded version of the above entangled measurements, we can per- 
form a SUM gate between any pair of encoded qudits anywhere in the computer. 
Note that the ancilla qudit can itself be anywhere in the computer; it need not 
be in the same block as either data qudit, or in the corresponding place in a 
different block. 

Given the SUM gate and the results of the previous section, we can perform 
the full Clifford group on the encoded data for any stabilizer code. This completes 
the proof up to step || This part of the proof is a significant improvement on the 
method used in @ . In that paper, it was necessary to introduce full ancilla blocks 
to perform a CNOT. Here, we need only a single logical ancilla qudit. In the 
case where a block may encode many qudits, this can be a major improvement. 
The price is that we must potentially perform entangled measurements on more 
than one block. This means we will have to use larger ancilla states for the 
measurement; this results in a greater potential for error, so we will have to 
repeat the measurement more times, and perhaps perform error correction a bit 
more often. However, in many situations, the total number of physical gates we 
need will decrease. 

Note that this procedure works just as well if the two logical qudits in- 
volved in the SUM are in blocks made up of different numbers of physical qudits. 
This means we can interact qudits encoded with different sorts of codes fault- 
tolerantly, or change the encoding of a single qudit without losing the protection 
against errors at any time. 



6 Completing the Universal Set of Gates 



The set of universal gates can be completed by adding the higher-dimensional 
analog of the Toffoli gate jl2| 

\a)\b)\c)^\a)\b)\c+ab) . (45) 

It turns out that a generalization of Shor's fault-tolerant construction of the 
Toffoli gate || will work here. 

Suppose we prepare a three-qudit ancilla in the +1 eigenstate of the three 
operators 

Mi = (X <g> I ® I) SUM(2 -> 3) , (46) 
M 2 = (I®X®T) SUM(1 -» 3) , (47) 
M 3 = (Kg Kg Z)PHASE(1,2)- 1 . (48) 

SUM(z — > j) is a SUM gate performed with the zth qudit as control and the jth 
qudit as target. PHASER, j) is the PHASE gate 

PHASE|a)|6) = u ab \a)\b) (49) 

performed on the ith and jth qudits. One important fact to note is that both of 
these gates are in the Clifford group. Since we have already constructed the Clif- 
ford group, this will enable us to also construct the Toffoli gate. The appropriate 
state is 

L4)=]>>}|6)| a &} . (50) 

a , 6 

Now, given three data qudits, perform inverse SUM gates (i.e., \a)\b) — > 
\a)\b — a)) from the first and second ancilla qudits to the first and second data 
qudits, respectively, and a SUM gate from the third data qudit to the third ancilla 
qudit. Now we measure the last three qudits, the original data qudits, in the bases 
Z, Z ', and X, respectively. After performing the appropriate corrections (which 
will likely involve gates from the Clifford group as well as the Pauli group), we 
are left with the data in the first three qudits, which were originally the ancilla 
qudits. It turns out that after these operations, a Toffoli gate has been performed 
on the data qudits. 

To construct the appropriate ancilla state \ A), we can again follow Shor. The 
various states 

a . b 

for j = 0, . . . , d — 1 are related to \A) by 

\Aj) = (Kg, Kg, X 1 ) \A) . (52) 
The states \Aj), like \A), are +1 eigenstates of Mi and M 2 , but M 3 \Aj) = uj j \Aj). 



Furthermore, note that 



= £|a>|6)|c) . (53) 

j — a,fe,c 

This last state is easily constructed as the Fourier transform of |000). Then by 
measuring M3 for this state, we can collapse the state into one of the states \ Aj). 
By applying the operator A _J , we get the state \A). 

To measure M 3 fault-tolerantly for a quantum code, we prepare the CAT 
state J2j \ jj ' ' ' j) (using the same number of qudits as in a block of the code). 
M3 is in the Clifford group, so we can perform it by some sequence of tranversal 
operations and measurements. By conditioning the appropriate operations for 
the ith qudit on the «th qudit of the CAT state, we can conditionally perform 
M3 on the code depending on the CAT state. Note that a conditional operation 
in the d-dimensional case means applying Af 3 when the control qudit is in the 
state \j). Once we have done this, when the code is in an eigenstate of M3 with 
eigenvalue uj s , the CAT state ends up in the state 

\CAT s ) = Y^u js \M---3} ■ (54) 
3 

The various states \CAT S ) are orthogonal to each other, and so can be distin- 
guished by an appropriate measurement. This therefore gives us a measurement 
of M3, completing the construction of the Toffoli gate and step || of the proof. 

Because everything we do is transversal, single qudit errors in the CAT state 
cannot become more than single qudit errors in any single block of the code. 
Naturally, after creating the CAT state we should verify it to make sure there 
are no correlated errors. In addition, a single qudit error in the CAT state could 
give us the wrong measurement result. Therefore, the measurement of M3 should 
be repeated in order to sufficiently increase our confidence in the result. 
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